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A  THREE-DIMENSIONAL  AXISYMMETRIC  CALCULATION  PROCEDURE  FOR  TURBULENT  FLOWS  IN  A  RADIAL  VANELESS  DIFFUSER 

Lawrence  F.  Schumann  ... 

Propulsion  Laboratory  *■’ 

AVSCOM  Research  and  Technology  Laboratories  , 

Lewis  Research  Center 
Cleveland,  Ohio 


ABSTRACT 

An  analytical  model  is  proposed  to  calculate  the 
three-dimensional  axisymmetric  turbulent  flowfield  in 
a  radial  vaneless  diffuser.  The  model  assumes  that 
the  radial  and  tangential  boundary  layer  profiles  can 
be  approximated  by  power  law  profiles.  Then,  using  the 
5  integrated  radial  and  tangential  momentum  and  continu- 

55  ity  equations  for  the  boundary  layer  and  corresponding 

inviscid  equations  for  the  core  flow,  there  results  six 
ordinary  differential  equation  in  six  unknowns  wnich 
can  easily  be  solved  using  a  Runge-Kutta  technique.  A 
model  is  also  proposed  for  fully  developed  flow.  The 
results  using  this  technique  have  been  compared  with 
the  results  from  a  three-dimensional  viscous,  axisym¬ 
metric  duct  code  and  with  experimental  data  and  good 
quantitative  agreement  was  obtained. 

NOMENCLATURE 

a  speed  of  sound 

A  area 

skin  friction  coefficient 
c  specific  heat  ratio  at  constant  pressure 

P  -  P, 


cpr  pressure  coefficient,  p - p — 


F  entrainment  function 

h  diffuser  half-width 

6  -  6 
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shape  factor. 
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H,  shape  factor,  — 
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exponent  in  radial  velocity  profile 
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PSTD 
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mass  flow 
Mach  number 

exponent  in  tangential  velocity  profile 
pressure 

standard  pressure,  101325  N/m 
radius 

r  -  r? 

dimensionless  radial  distance,  - =~ 

r3  ~  r2 

gas  constant 
Reynold's  number 
entropy 

streamline  distance 
temperature 

standard  temperature,  288.2  K 
velocity 

coordinate  normal  to  the  wall 
dimensionless  y-distance,  y/h 
flow  angle  from  radial  direction 
specific  heat  ratio 
boundary  layer  thickness 

,« 

displacement  thickness. 


,/  /l-^- 

J0  {  vre 


dy 


”  Ik  ' 

.  %m  m  SLm 


-«  .  ■  ^  J 


*.  *v  \  \  *,  •.  •,  \ 

■  i.m  Cm  ’Ck  K ■  V.  lT*  Urn l 


im  XmJL* 


«ls  displacement  thickness,  J  ^l--^-jdy 

momentum  thickness,  J"*  -f—  f  1  -  tr~ ^dy 
‘  '  0  vre  \  vre  / 


«2S  momentum  thickness,^  (l  --^-\dy 

u  coefficient  of  viscosity 

p  density 

t  shear  stress 

»  entropy  function,  J*  c^ 

Subscripts: 
c  centerline 

cr  critical 

e  edge  of  boundary  layer 

r  radial  direction 

R  reference  condition 

s  streamwise  direction 

w  wall 

y  direction  normal  to  the  wall 

9  tangential  direction 

1  compressor  inlet 

2  vaneless  diffuser  inlet 

3  vaneless  diffuser  exit 
Superscripts: 

'  total  conditions 

(  )  mass-averaged  conditions 

INTRODUCTION 

The  flowfield  in  the  radial  vaneless  diffuser  of  a 
centrifugal  compressor  is  extremely  complex  since  the 
flow  is  turbulent,  unsteady,  viscous,  and  three- 
dimensional.  Also,  depending  on  the  initial  state  of 
the  end-wall  boundary  layers  and  the  diffuser  length, 
the  flow  may  become  fully  developed  (i.e.,  the  end-wall 
boundary  layers  may  merge  in  the  center  of  the  channel) 
or  may  separate  off  one  of  the  walls.  One  of  the  ear¬ 
liest  attempts  to  calculate  this  flowfield  was  due  to 
Stanitz  (1).  Other  investigators  have  used  a  similar 
approach  to  calculate  the  frictional  losses  in  the 
vaneless  diffuser  (2,3).  In  Stanitz'  model,  the  flow 
is  assumed  to  be  one-dimensional  with  the  losses 
accounted  for  by  using  a  shear  force  term  proportional 
to  an  assigned  friction  factor  in  the  radial  and  tan¬ 
gential  momentum  equations.  This  approach  has  several 


weaknesses.  First  of  all,  the  development  of  the  end- 
wall  boundary  layers  is  not  accounted  for  except 
through  the  inclusion  of  an  effective  channel  height 
term.  This  generally  leads  to  an  overprediction  of  the 
static  pressure  rise  unless  an  unrealistically  high 
friction  factor  is  used.  Secondly,  the  method  gives 
no  information  about  the  velocity  profiles  and  flow 
angle  distribution  at  the  exit  of  the  vaneless  dif¬ 
fuser.  This  information  would  be  helpful  in  the  design 
of  the  leading  edge  of  the  diffuser  vanes.  Finally, 
the  choice  of  an  appropriate  friction  factor  and  the 
effective  passage  height  distribution  are  totally 
arbitrary  and  will  often  vary  from  one  machine  to 
another. 

More  recently,  several  investigators  (4-6)  have 
attempted  to  calculate  the  three-dimensionallioundary 
layers  in  vaneless  diffusers.  The  major  drawback  of 
these  methods  is  the  empiricism  used  in  the  specifica¬ 
tion  of  the  velocity  profiles.  This  affects  the  cal¬ 
culation  of  the  boundary  layer  shape  factor  and 
therefore  also  the  skin  friction  and  the  location  of 
possible  separation  points. 

Currently,  there  are  computer  codes  (7)  which  can 
calculate  the  three-dimensional  flow  fields  in  annular 
ducts.  These  codes  can  also  be  used  for  calculating 
the  flow  field  in  a  radial  vaneless  diffuser  but  their 
running  times  and  complexity  preclude  their  use  as 
design  tools.  Also,  in  the  case  of  the  ADD  code  (7) , 
there  are  stringent  Mach  number  and  Reynold's  number 
1  imitations. 

In  the  present  paper,  a  three-dimensional  axisym- 
metric  calculation  procedure  for  turbulent  flows  in  a 
radial  vaneless  diffuser  is  presented.  The  method  is 
a  noniterative  integral  method  and  is  therefore  much 
faster  than  finite  difference  methods.  The  primary 
assumption  made  is  that  both  the  axial  and  tangential 
boundary  layer  velocity  profiles  can  be  described  by 
power  laws;  however,  the  exponents  in  these  power  laws 
are  allowed  to  vary  with  radius  as  prescribed  by  the 
equations  of  motion.  Compressibility  is  accounted  for 
in  the  radial  direction  but  not  in  the  direction  normal 
to  the  end-walls.  Turbulent  flow  is  assumed  everywhere 
using  the  semi-empirical  Ludwieg-Ti liman  relation  for 
skin-friction.  A  method  for  calculating  fully- 
developed  flows  is  also  presented.  Comparisons  are 
made  to  results  from  the  ADD  code  for  both  the  boundary 
layer  analysis  and  for  the  fully  developed  flow  analy¬ 
sis.  Comparisons  are  also  made  with  experimental  data 
in  the  literature.  In  all  cases,  good  agreement  is 
demonstrated. 

ANALYSIS 

The  type  of  vaneless  diffuser  geometry  considered 
by  this  method  of  analysis  is  shown  in  Fig.  1.  Typical 
velocity  profiles  for  the  cases  of  boundary  layer  flow 
and  fully  developed  flow  are  shown  in  Figs.  2(a)  and 

Boundary  Layer  Flow 

The  following  assumptions  are  made  for  this  analy¬ 
sis: 

(1)  The  flow  is  steady  and  turbulent  everywhere. 

(2)  The  fluid  is  a  perfect  gas. 

(3)  The  flow  is  axi symmetric. 

(4)  The  end-walls  of  the  vaneless  diffuser  are 
adiabatic. 

(5)  There  is  no  density  variation  in  the  direction 
normal  to  the  vaneless  diffuser  end-walls  although  the 
density  is  allowed  to  vary  in  the  radial  direction, 
i.e. ,  p  =  p(r) . 

(6)  The  vaneless  diffuser  geometry  and  flow  are 
symmetric  about  the  diffuser  centerline. 


(7)  The  boundary  layer  velocity  profiles  in  the  dv  v 

radial  and  tangential  directions  may  be  described  by  — r- ^  =  — —  (10) 

the  power  law  equations  ar  r 


where  vrp,  v9e,  a,  m,  and  n  are  functions  of 
radius  only. 

The  boundary  layer  equations  of  motion  in  cylin¬ 
drical  coordinates  for  axisymmetric  turbulent  flow  are 
as  follows: 


From  the  boundary  layer  assumption,  which  is  inherent 
in  Eq.  (5),  the  static  pressure  is  constant  across  the 
channel  and  thus  is  a  function  of  radius  only. 

The  remaining  three  equations  are  obtained  from 
channel  continuity,  the  energy  equation,  and  the  equa¬ 
tion  of  state.  Channel  continuity  can  be  written 

">  =  2pvre(2»r)(h  -  «lr)  (11) 


Continuity 


where  for  a  power  law  profile 
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(2)  5lr  =  r^i< 

Differentiating  Eq.  (11)  with  respect  to  the  radius 
yields  an  expression  for  dvre/dr.  If  the  total  tem¬ 
perature  is  assumed  to  be  constant  in  the  vaneless 
diffuser,  which  is  true  if  there  are  no  inlet  normal 
/ 3 \  gradients  of  total  temperature,  the  energy  equation  for 
'  the  core  flow  can  be  written 


Tangential  Momentum 
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(5) 


If  the  power  law  profiles  of  Eq.  (1)  are  assumed, 
the  problem  reduces  to  that  of  determining  the  free 
stream  velocity  components  vre  and  vee;  the 
boundary  layer  thickness,  a;  the  boundary  layer  profile 
exponents,  m  and  n;  and  the  free  stream  thermodynamic 
properties,  p,  p,  and  Te.  Thus  eight  equations  are 
required  to  completely  solve  the  flow  field.  Three  of 
these  equations  are  obtained  by  integrating  Eqs.  (2)  to 
(4)  through  the  boundary  layer  which  results  in  the 
fol lowing 


The  final  equation  which  completes  the  system  is  the 
equation  of  state  for  the  core  flow  which  is 

p  =  pRgre  (13) 

One  additional  equation  is  needed  to  solve  this 
system  since  an  additional  unknown  is  added  in  the 
integration  of  Eqs.  (2)  to  (4);  i.e.,  the  normal  veloc¬ 
ity  at  the  edge  of  the  boundary  layer,  vye.  This  can 
be  modelled  using  Head's  entrainment  relation  which 
empirically  predicts  the  rate  at  which  fluid  from  the 
free  stream  is  entrained  into  the  boundary  layer.  Head 
showed  that  the  entrainment  rate  can  be  modelled  as  a 
function  only  of  a  boundary  layer  shape  factor.  Sumner 
and  Shanebrook  (8)  extended  the  definition  to  three- 
dimensional  compressible  flow  as  follows 
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Equations  (3,4)  can  also  be  evaluated  for  the 
isentropic  core  flow.  This  results  in  the  following 
two  equations 


F=f  (veS?-  Vye)  =  °.°306  (Hi  -  3) 


-0.653 


(14) 


Finally,  expressions  for  the  wall  shearing  stresses 
can  be  given  by 


pVre  cf 

Trw  cos  a_  ~1 
e 


pvre  .  cf 

Tew  =  cosIT  tan  “e  1 
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(15) 


(16) 


where  Cf  is  given  by  the  Ludwieg-Tillman  relation 
which  is 

c f  =  0.246  r^0-268  ^  -1'561  H2  (17) 
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(9) 


Using  the  velocity  profiles  from  Eq.  (1)  and  per¬ 
forming  the  indicated  differentiations  and  integra¬ 
tions,  Eqs.  (6)  to  (8)  and  (11)  can  be  combined  to 
yield  the  following  equations 
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+  (20) 
l  +  m  Ur  cos  n  '  1 


dn  _  n  dm  + 

Ur  =  1  +  m  clr 


(1  +  m  +  n)2  T  Cf  nF  1 

7'cosTT -  L“?“T1  +'n>  +  n)  J  (21> 


Equations  (9),  (10),  and  (18)  to  (21)  can  then  be 
solved  sequentially  for  the  six  unknowns  (vre,  m,  t,  n, 
p,  and  v„e)  using  a  fourth  order  Runge-Kutta  tech¬ 
nique  sucn  as  that  described  in  Ref.  (9).  Equation 
(18)  is  good  only  for  symmetric  vaneless  diffusers  with 
moderate  wall  curvatures  and  uniform  core  flows.  How¬ 
ever,  Eqs.  (9),  (10),  and  (19)  to  (21)  are  good  for 
nonuniform  core  conditions  and  could  be  solved  together 
with  the  equations  describing  the  core  flow  as  an 
interactive  boundary  layer  problem.  The  above  system 
of  equations  is  aiso  good  for  incompressible  flow  if 
the  Mach  numbers,  Mre  and  Mge,  are  set  equal  to 
zero  in  Eqs.  (18)  and  (20). 

A  complete  derivation  of  the  final  form  of  the  con¬ 
tinuity  equation  is  described  in  the  Appendix.  The 
radial  and  tangential  momentum  equations  are  derived 
simi larly. 


JL.  Jl')  •  V  - 
vrc  'h'  ’  V  ^ 


where  vr(~,  ve(-,  h,  m,  and  n  are  functions  of  radius 
only.  (3)  Analagous  to  fully  developed,  turbulent, 
pipe  flow  theory  (10),  the  shear  stress  is  assumed  to 
vary  linearly  across  the  channel. 

For  fully  developed  flow,  an  isentropic  core  flow 
no  longer  exists.  Therefore  Eqs.  (9)  and  (10)  are  no 
longer  valid.  However,  with  the  assumption  of  a  linear 
variation  of  shear  stress  across  the  channel,  Eqs.  (3) 
and  (4)  can  be  evaluated  along  the  channel  centerline 
which  yields  the  equations 


rc  . 

p  V  — 3 + 

u  rc  dr 


vec  _  _1 _ ’ew 

r  -  Pvrc  "FT 


where  xrw  and  are  given  by  Eqs.  (15)  and  (16) 

and  Cf  is  given  Dy  Eq.  (17)  with 


h2  -  1 


If  the  continuity  equation  is  integrated  across  the 
vaneless  diffuser  half-width,  the  following  equation 
results 


|  .r,r 


dy  =  0 


For  the  power  law  profiles  of  Eq.  (22),  this  becomes 
d  v 

1  ^  1  dp  A  1  rc  ,  1  dh  1  dm  n 

The  radial  momentum  equation,  Eq.  (3),  can  also  be 
integrated  across  the  diffuser  half-width.  This  be¬ 
comes,  after  using  the  results  from  Eqs.  (23)  and  (26) 
and  the  velocity  profiles  from  Eq.  (22) 


37  =  -  2m(l  ♦  m)  ( 1  ♦  2m) 


dvrc  n(  1  ♦  2m)  tanV| 


When  the  tangential  momentum  equation  is  integrated 
across  the  diffuser  half-width,  the  following  equation 
results 


Fully  Developed  Flow 

The  system  of  equations  described  in  the  previous 
section  for  boundary  layer  flow  breaks  down  when  the 
boundary  layers  on  the  two  end-walls  meet  in  the  center 
of  the  channel,  i.e.,  when  «  j>  h.  In  this  case,  the 
flow  is  said  to  be  fully  developed  and  the  velocity 
profiles  are  assumed  to  be  as  shown  in  Fig.  2(b). 

The  following  assumptions  are  made  for  fully  devel¬ 
oped  flow  in  addition  to  those  already  made  for  bound¬ 
ary  layer  flow.  (1)  The  boundary  layer  equations, 

Eqs.  (2)  to  (5),  remain  valid  for  fully  developed  flow. 
(2)  The  velocity  profiles  can  be  described  by  the  power 
law  equations 


dn 

dr 


_n _ dm  +  (m  +  n)(l  +  m  +  n)  f 

+  m  dr  cos  a,  Jh 


(28) 


The  density  variation  in  the  throughflow  direction 
can  be  obtained  from  a  combination  of  the  equation  of 
state  and  the  energy  equation.  Then,  using  the  above 
relations,  the  following  equation  for  the  variation  of 
the  centerline  radial  velocity  component  results. 
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2h  COS  a. 


c 

(29) 


Equation  (29)  together  with  Eqs.  (23),  (24),  (27), 
and  (28)  then  forms  a  system  of  five  ordinary  differen¬ 
tial  equations  in  five  unknowns  (vrc,  v9c,  m,  n, 
and  p)  which  can  be  solved  using  a  fourth-order  Runge- 
Kutta  technique. 


Separated  flow 

By  examining  Eqs.  (19)  and  (27),  it  can  be  seen 
that  the  radial  velocity  profile  exponent,  m,  and 
therefore  also  the  shape  factor,  H2,  will  increase 
rapidly  if  there  is  a  large  rate  of  deceleration  or  if 
the  flow  angle  becomes  large.  This  usually  means  that 
flow  separation  is  imminent.  A  large  deceleration  rate 
is  caused  by  a  rapid  increase  in  flow  area.  A  large 
flow  angle  can  be  the  result  of  the  impeller  design  but 
will  more  often  be  the  result  of  a  low  mass  flow  rate 
through  the  compressor.  In  the  latter  case,  a  large 
flow  separation  may  occur  in  the  vaneless  space  which 
could  precipitate  compressor  stall  and  subsequent  surge 
of  the  system. 

The  methods  outlined  in  the  previous  two  sections 
for  boundary  layer  flow,  and  fully  developed  flow, 
could  be  used  to  calculate  through  a  separated  region 
of  the  flow.  However,  if  this  were  done,  large  values 
of  the  velocity  profile  exponent,  m,  would  be  calcu¬ 
lated  which  would  result  in  unrealistic  velocity  pro¬ 
files.  A  possible  way  of  treating  this  problem  would 
be  to  assume  a  constant  value  for  the  exponent,  m, 
after  separation.  Separation  could  be  determined  by 
examining  the  shape  factor,  Hj,  or  the  skin  friction 
coefficient,  Cf.  The  problem  could  then  be  treated 
in  the  inverse  mode  using  either  Eq.  (19)  for  boundary 
layer  flow  or  Eq.  (27)  for  fully  developed  flow  to 
determine  the  variation  of  free-stream  radial  velocity. 
However,  since  power  law  profiles  cannot  describe  neg¬ 
ative  velocities  which  occur  when  the  flow  separates, 
it  would  probably  be  necessary  to  use  some  other  formu¬ 
lation  for  the  separated  flow  problem. 


meridional  velocity  component  in  the  impeller.  How¬ 
ever,  s  will  probably  be  a  significant  portion  of  the 
impeller  exit  half-width. 

If  the  compressor  mass  flow  is  specified,  the  core 
flow  radial  velocity  component  can  then  be  calculated 
from  Eq.  (11).  However,  the  core  flow  tangential 
velocity  component  and  corresponding  velocity  profile 
exponent  are  more  difficult  to  determine.  A  boundary 
layer  analysis  of  the  impeller  would  give  values  for 
the  relative  tangential  velocity  component  and  its 
corresponding  exponent.  However,  as  the  fluid  moves 
from  the  relative  frame  of  reference  of  the  impeller 
(assuming  a  shrouded  impeller)  to  the  absolute  frame 
of  the  vaneless  diffuser,  a  large  shear  occurs  on  the 
fluid  as  it  attempts  to  decelerate  from  a  velocity  near 
the  wheel  speed  to  zero  near  the  end-walls.  Thus  it  is 
impossible  to  model  the  velocity  profile  in  this  region 
by  a  power  law.  Experimental  data  does  show,  however, 
that  a  velocity  profile  that  can  be  approximated  by  a 
power  law  is  quickly  established.  Experimental  data 
also  shows  (4j  that  an  exponent  of  n  0.1  adequately 
describes  the  tangential  boundary  layer  at  the  vaneless 
diffuser  inlet.  The  core  flow  tangential  velocity  com¬ 
ponent  can  then  be  obtained  from  a  conservation  of  tan¬ 
gential  momentum  analysis.  In  any  event,  the  initial 
values  of  the  exponents,  m  and  n,  do  not  have  a  large 
effect  on  the  final  solution. 

For  the  case  of  fully  developed  flow  entering  the 
vaneless  diffuser,  the  same  variables  must  be  specified 
as  for  the  boundary  layer  case  with  the  exception  of 
the  boundary  layer  thickness  (which  is  equal  to  the 
diffuser  half-width,  h).  For  this  case,  the  value  of 
as  in  Eq.  (30)  should  also  include  the  additional  core 
flow  loss  due  to  a  nonzero  shear  stress  gradient  at  the 
centerline.  The  values  of  the  other  variables  are  then 
calculated  in  the  same  way  as  for  boundary  layer  flow. 

Loss  Calculation 

The  loss  in  a  vaneless  diffuser  can  be  expressed 
most  easily  and  most  rigorously  by  an  increase  in  the 
mass-averaged  entropy  of  the  fluid.  The  mass-averaged 
entropy  at  a  specified  radius  in  the  vaneless  diffuser 
can  be  expressed 


_  2p  JQ  (s  -  s?e)  vr(2tr)dy 

(s '  '*) - l 


Initial  Conditions 

In  order  to  solve  Eqs.  (9),  (10),  and  (18)  to  (21) 
for  boundary  layer  flow,  initial  values  for  the  six 
variables  (p,  i,  vre,  v9e,  m,  and  n)  must  be 
specified.  The  value  of  the  static  pressure  can  be 
obtained  from  experimental  data  or  from  an  analysis  of 
the  impeller  flowfield.  In  the  latter  case,  the  pres¬ 
sure  could  be  obtained  from  the  equation 


where  the  local  entropy  is  given  by 

T 


—  j  —  AS 

p  =  P1  e  *g  (30) 

The  entropy  rise,  as,  in  Eq.  (30)  should  include  the 
entropy  rise  due  to  the  mixing  of  the  impeller  blade 
surface,  boundary  layers,  wakes,  and  any  other  losses 
which  are  assumed  to  affect  the  impeller  core  flow. 

The  values  of  «  and  m  can  only  be  obtained  from 
an  analysis  of  the  impeller  end-wall  boundary  layers. 
Ordinarily  the  value  of  m  should  be  on  the  order  of 
0.2  or  0.3  if  there  is  no  great  deceleration  of  the 


The  entropy  rise  through  the  diffuser  can  then  be 
expressed 


(s3  -  s2)=  (s  -  s2e)3  -  (s  -  s2e)?  (33) 

where  the  values  of  (s  -  S2?)2  and  (s  -  s2e)3 
are  obtained  by  a  numerical  integration  of  Eq.  (32) 
across  the  boundary  layer.  Since  the  total  temperature 
in  the  vaneless  diffuser  is  assumed  to  be  constant 
everywhere,  the  mass-averaged  total  pressure  ratio  can 
be  expressed 
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ordinary  differential  equations.  The  continuity  equa¬ 
tion  will  serve  as  an  example  of  how  this  is  done. 

Also,  the  continuity  equation  shows  how  the  entrainment 
function  is  used  in  the  derivation.  If  Leibnitz'  rule 
is  used  to  interchange  differentiation  and  integration, 
Eq.  (6)  becomes 


The  term  on  the  right  hand  side  of  Eq.  (35)  is  Head's 
entrainment  function,  F,  and  will  be  evaluated  later. 
Since  the  density  is  assumed  to  be  a  function  of  radius 
only,  it  comes  out  of  the  integral  and  the  integrals 
can  be  evaluated  for  the  power  law  profiles  of  Eq.  (1) 
which  results  in  the  following  equation 

1  d«  +  _1_  dvre  +  1  dp  +  \  _  1  dm  _  F(1  +  m) 

«  Ur  vre  dr"  p  Ur  r  1  +  m  Ur  ~  cos  ag 

(36) 

This  is  the  form  desired  except  for  the  density  term. 
The  density  variation  can  be  expressed  in  terms  of  the 
pressure  variation  for  the  isentropic  core  flow  as 
follows 


dp  1  dp 

Ur  =  "7  dr 

d 

Combining  this  with  Eq.  (9)  then  yields 


(37) 
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(39) 


Using  similar  techniques,  the  radial  and  tangential 
momentum  equations  can  also  be  reduced  to  ordinary 
differential  equations. 

The  entrainment  function,  F,  is  given  by  Eq.  (14). 
The  shape  factor.  Hi,  is  defined  by  the  expression 


H 


1  = 


(40) 


De  Ruyck  and  Hirsch  (11)  approximated  this  using  the 
radial  velocity  components  in  the  definitions  of  «i 
and  «2»  thus 


H 


1  - 


(41) 


Although  Deruyck  and  Hirsch  stated  that  Eq.  (41)  is 
good  only  for  small  angles,  a,  Davis  ( 6 )  shows  that  the 
radial  thicknesses,  «ir  and  6£r,  are  the  proper 
thicknesses  to  use  in  the  definition  of  Hi  for  use 
in  the  entrainment  function.  Using  the  power  law  pro¬ 
files  of  Eq.  (1)  in  the  definition  of  «ir  and  «2r» 
the  shape  factor.  Hi,  then  becomes 


1_  dp^ 

p  dr 


M  dv 
re  re 

vre 


(38) 


1  +  2m 
1  m 


(42) 


and  finally  the  entrainment  function,  F,  becomes 


The  final  form  of  the  continuity  equation  then  becomes 


0.653 


(43) 
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TABLE  I.  -  FIVE  TEST  CASES  FOR  COMPARISON  OF  ANALYSIS 

TO  RESULTS  FROM  THE  ADD  CODE 

[Parallel  wall  diffusers  -  r2  =  8.049  cm,  r2/r2  = 
1.1605,  h/r2  =  0.03329,  T'/TST0  =  1.] 


Inlet  conditions 


Case 

6/h 

m 

n 

P/PSTD 

vre/vcr 

vec/vcr 

“e 

0.2400 

0.2 

0.2 

0.7087 

0.7498 

0 

0 

.2400 

.2 

.2 

.9172 

.2705 

0.2705 

Em 

.9556 

.2 

.2 

.9172 

.2705 

.2705 

m 

.2400 

.2 

.2 

.7009 

.2603 

.7151 

70 

.9556 

.2 

.2 

.7009 

.2603 

.7151 

70 

TABLE  II.  -  TWO  TEST  CASES  FOR  COMPARISON  WITH  DATA  FROM 

REFERENCE  4 


[Parallel  wall  diffusers  -  r2  =  17.64  cm,  r2/r2  = 
1.377,  h/r2  =  0.0605,  T7TSTD  =  !•] 


Inlet 

conditions 

Case 

6/h 

m 

n 

P/PSTD 

Ve» 

m/sec 

Vee» 

m/sec 

“e 

VI 

0.288 

m 

HI 

1 

6.83 

1 

46.5 

VII 

.360 

\m 

MM 

1 

6.35 

67.4 

PRESENT  ANALYSIS 
ADD  CODE 


.2  .4  .6  .8  1. 

Y 

(a)  Case  I:  Boundary  layer  flow*,  02e  =  0°. 

Figure  3.  -  Velocity  profile  comparisons  with 
results  from  ADD  code. 
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-  PRESENT  ANALYSIS 

- ADD  CODE 

O  NUMERICAL  INSTABILITY 


(0  Case  IV.  (d)  Case  V. 

Figure  5.  -  Centerline  flow  angle  comparison  with 
ADD  code. 
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